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' Abstract. We present maps of the 850 fim and 450 fim continuum emission seen towards a sample of 68 high-mass 

3^ , protostellar candidates with luminosities ranging from lO'^'^ Lq to ~ 10^ Lq . Most of these candidate high-mass 

stars are in the earliest stages of evolution, and have not yet developed an ultra-compact Ifll region. We observe a 
1/^ ' variety of continuum emission morphologies, from compact symmetric sources through to multiple cores embedded 

in long filaments of emission. We find on average there is a 65% probability of an IRAS point-source having a 
companion detection at submillimetre wavelengths. The ratio of integrated flux to peak flux for our detections 
^ I shows no strong dependence on distance, suggesting the emission we have observed is primarily from scale-free 

. envelopes with power-law density structures. Assuming a near kinematic distance projection, the clumps we detect 

vary in mass from ~1 M0 to over 1000 M0, with a mean clump mass of 330 Mq, column density of 9 x 10'^'' 
cm~^ and diameter of ~ 0.6 pc. The high luminosity and low mass of the smallest clumps suggests they are 
^ ■ accompanied by a minimal number of stellar companions, while the most massive clumps may be examples of 

' young protogroups and protoclusters. We measure the spectral index of the dust emission (a) and the spectral 

, index of the dust grain opacity (/3) towards each object, finding clumps with morphologies suggestive of strong 

temperature gradients, and of grain growth in their dense inner regions. We find a mean value for f3 of 0.9, 
_ significantly smaller than observed towards UCHII regions. 
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I years, and the concepts of gravity-driven collapse and 
■ - - ' accretion-driven evolution appear to consistently explain Unfortunately, a powerful protostar and UCHII region 

how a low-mass pre-main-sequence star can form from soon act to disrupt and confuse their surroundings, so the 

a cloud core. However, our knowledge of how high-mass initial conditions of the natal cloud and the mechanisms 

stars form has remained limited, primarily due to a lack that led to the formation of the massive protostar can- 

of candidate high-mass protostars to study. At the up- not be unambiguously reconstructed. As a result, many 

per reaches of the initial mass function, high-mass stars questions about high-mass protostars remain - in partic- 

are statistically rare, and coupled with a characteristically nlar, do they form via processes similar to their low-mass 

short evolutionary timescale, there have been few chances counterparts? To address the mechanisms that create and 

to observe massive stars at the instance of formation. shape high-mass stars, we must observe before they have 

Until recently, most young high-mass stars were first formed a UCHII region, during the initial collapse of the 

identified through the detection of a radio-bright ultra- star-forming core, 
compact HII (UCHII) region, considered a beacon point- 
ing to the presence of a young high-mass star. As high- 
mass protostars increase in mass and luminosity, they emit 

an ever larger number of high energy UV photons which ^, 1 • ii ^ i ■ 1 ^ j_ 1 

° °, ,. , Clumps bearing the youngest high-mass protostars have 
ionize the protostars immediate surroundings, hence the 1 j.- 1 1 j-jn ^ c j 1 • 
° proved particularly ditncuit to hnd, even though their 

Send offprint requests to: G. A. Fuller, e-mail: identifying characteristics have been known for almost 25 

g.fuller@umist.ac.uk years. For example, Habing & Israel (1979) predicted that 



1. Introduction small, compact nature of a UCHII region is usually consid- 

1 ■ r 1 • 1 1 • 1 eied evidence of the youthful status of the driving source 

Our understanding 01 the processes involved m low-mass / ,,1 u j u j- j.- u ^ j. 1 r 

_ . , ^, , 1 n (although debate continues about the exact timescale 01 

star formation has matured steadily over the last forty ^^r^^r^^ ^ r ■ ^ x 1 nnnn\ 

, , p . , . , the UCHII stage; for a review, see Kurtz et al. 2000). 



1.1. The search for precursors of UCHII regions 
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candidate high-mass protostars should be founded embed- 
ded in dense environments, and although highly luminous, 
they should not at this stage be associated with HII re- 
gions. However, despite knowledge of these distinguishing 
features, it was not until the last decade that samples of 
candidate high-mass protostellar objects (HMPOs) were 
finally compiled. 

Observations of these preliminary samples of HMPOs 
have allowed the first glimpses of high-mass protostars in 
their earliest evolutionary states: a typical core not yet as- 
sociated with an ionised region is found to be larger, more 
massive, and more turbulent than a UCHII-class proto- 
star, with a typical diameter of around 0.5-1.0 pc and a 
mass that may range from a few tens to a few thousand 
solar masses (Brand et al. 2001; Beuther et al. 2002a). 
They are cooler, with typical dust temperature averag- 
ing around 30-40 K (Sridharan et al. 2002; Molinari et al. 
2000), while the dust opacity usually has a spectral index 
of around 2, suggestive of silicate dust grains (Molinari et 
al. 2000). Self-absorption profiles towards a number of can- 
didate HMPOs suggest infall may be an important part of 
the formation mechanism (eg. Brand et al. 2001; Fuller et 
al. 2004), while outfiow observations suggest that accre- 
tion is a significant process (Zhang et al. 2001; Beuther et 
al. 2002b; Molinari et al. 2002). Water maser emission has 
also been detected towards candidate protostars, a feature 
thought to be missing from more evolved sources (eg. Palla 
et al. 1993; Sridharan et al. 2002). Despite these advances, 
there is still much to be learned about the pre-UCHH stage 
of high-mass star formation, and there remains a need for 
additional candidates and further observations. 

1.1.1. A new sample of high-mass protostars 

Recently, Sridharan et al. (2002; SBSMW hereinafter) 
identified a new sample of HMPOs. The SBSMW sample 
is a flux-limited sample, constructed through an analysis 
of the IRAS point-source catalogue: as young high-mass 
stars are usually associated with UCHH regions, they be- 
gan by initially selecting bright IRAS detections {Sqq > 90 
Jy and Sioo > 500 Jy) with colour characteristics similar 
to known UCHH regions (they conform to the Wood & 
ChurchweU (1989) FIR colour criteria that selects UCHII 
regions, and they also satisfy the additional Ramesh & 
Sridharan (1997) criteria). Candidate sources detected 
in Galaxy-wide 5GHz continuum surveys were removed, 
thereby rejecting sources already sufficiently evolved to 
have ionized their surroundings. As a final requirement, 
successful candidates must also be associated with CS(2-1) 
emission, an indicator of dense molecular gas (Bronfman 
et al. 1996). 

In total, sixty-nine IRAS point sources satisfied these 
cumulative criteria, identiiying these sources as poten- 
tially among the most massive and deeply embedded pre- 
UCHII protostars in our Galaxy. The SBSMW sample has 
been studied in detail over the last few years, and their 
status as high-mass candidate protostars has been sup- 



ported through observations of 1.2mm and 3.6cm contin- 
uum emission (Beuther et al. 2002a; SBSMW), molecular 
line emission (CS, CO and NH3), and H2O and CH3OH 
maser transitions (SBSMW) towards the sources. 

This paper presents the results of a new set of sub- 
millimetre (submm) observations of the SBSMW sam- 
ple of candidate high-mass protostars. All but one (IRAS 
18517-F0437) of the SBSMW candidate HMPOs were ob- 
served. An additional source. IRAS 18449-0158, was ob- 
served but this source does not satisfy the SBSMW criteria 
and is not included in any analysis. Our observations are 
detailed in §2, with maps of the reduced data foimd in §3. 
We measure the multiplicity of the detections in §3.1, com- 
menting on the position and morphologies of the sample 
in §3.4. We analyse the dust optical depth in §4, and use 
the spectral index of the emission to investigate the nature 
of the dust in §4.1. We calculate the mass characteristics 
of our sample in §5. and consider the implications of the 
cumulative mass spectrum in §5.1. After a brief discussion 
and comparison of our results with the IRAM 1.2mm con- 
timmm observations of Beuther et al. (2002a) in §6.1, we 
conclude in §7 with a summary of our results. 

This paper presents the first half of our study and anal- 
ysis of the dust emission; the companion to this paper 
presents the results of radiative transfer modelling of the 
clumps (Williams et al. 2004). 

2. Observations and data reduction 

The sample of HMPOs was observed at 850 /im and 
450 /im between March 2000 and June 2000 using the 
Submillimetre Common-User Bolometer Array (SCUBA) 
on the James Clerk Maxwell Telescope (Holland et 
al. 1999). The SBSMW sample target co-ordinates and 
date(s) of observation are listed in Table 1. The SCUBA 
array covers a hexagonal 2.5' field of view with 97 and 
37 pixels at 450 fim and 850 /im respectively. Maps were 
formed simultaneously at both frequencies using the "jig- 
gle" mode, in which the telescope beam is moved around a 
64-point pattern by the secondary mirror in order to fully 
sample the sky. 

The data were reduced using the SCUBA User 
Reduction Facility (SURF; Jenness & Lightfoot 1998). 
Correlated sky noise was removed using the REMSKY 
routine, based on the signal from a hand-picked sam- 
ple of bolometers considered free from source emission. 
Maps were extinction calibrated from skydips and flux 
calibrated in terms of Jy beam~^ from maps of Uranus, 
IRAS 16293-2422, CRL 618, and CRL 2688, following the 
procedures defined by Sandell et al. (2001). 

Zenith opacities at 225 GHz ranged from 0.05-0.12 dur- 
ing the observations, but usually averaged around 0.10. 
Telescope pointing was calibrated many times during each 
observing run, and telescope drift was minimal, requiring 
very small {a = 1.7") corrections overall. We measured 
the JCMT beam size from observations of Uranus, find- 
ing a full-width half-maximum of 9beam = 8.0" at 450 /xm 
and 6beam = 14.4" at 850 fim. An average 1-a RMS noise 
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IRAS Source 



Position (J2000) 
a S 



Distance (kpc) 
far near 



Date(s) of 
observation 



05358+3543 


05 


39 


10, 


.4 


+35 


45 


19 


1 


.8 




03/18 


05490+2658 


05 


52 


12, 


,9 


+26 


59 


33 


2 


.1 




03/18 


05553+1631 


05 


58 


13, 


,9 


+ 16 


32 


00 


2 


.5 




03/18 


18089-1732 


18 


11 


51, 


,3 


-17 


31 


29 


13 


.0 


3.6 


05/11 


18090-1832 


18 


12 


01, 


,9 


-18 


31 


56 


10 


.0 


6.6 


05/11 


18102-1800 


18 


13 


12, 


,2 


-17 


59 


35 


14 


.0 


2.6 


05/11 


18151-1208 


18 


17 


57, 


,1 


-12 


07 


22 


3 


.0 




05/11 


18159-1550 


18 


18 


47, 


.6 


-15 


48 


54 


11 


.7 


4.7 


05/11 


18182-1433 


18 


21 


07, 


,9 


-14 


31 


53 


11 


.8 


4.5 


05/11 


18223-1243 


18 


25 


11, 


,1 


-12 


42 


15 


12 


.4 


3.7 


05/11 


18247-1147 


18 


27 


31, 


,1 


-11 


45 


56 


9 


.3 


6.7 


05/11 


18264-1152 


18 


29 


14, 


,3 


-11 


50 


26 


12 


.5 


3.5 


05/11 


18272-1217 


18 


30 


02, 


7 


-12 


15 


27 


2 


.9 




05/11 


18290-0924 


18 


31 


44, 


,8 


-09 


22 


09 


10 


.5 


5.3 


05/11 


18306-0835 


18 


33 


21, 


,8 


-08 


33 


39 


10 


.7 


4.9 


05/11 


18308-0841 


18 


33 


31, 


,9 


-08 


39 


17 


10 


.7 


4.9 


05/11 


18310-0825 


18 


33 


47, 


.2 


-08 


23 


35 


10 


.4 


5.2 


05/11 


18337-0743 


18 


36 


29, 


,0 


-07 


40 


33 


11 


.5 


4 


05/11, 


18345-0641 


18 


37 


16, 


,8 


-06 


38 


32 


9 


.5 




05/11 


18348-0616 


18 


37 


29, 


,0 


-06 


14 


15 


9 


.0 


6.3 


05/11 


18372-0541 


18 


39 


56, 


,0 


-05 


38 


49 


13 


.4 


1.8 


05/11 


18385-0512 


18 


41 


12, 


,0 


-05 


09 


07 


13 


.1 


2 


05/23 


18426-0204 


18 


45 


12, 


,8 


-02 


01 


12 


13 


.5 


1.1 


05/23 


18431-0312 


18 


45 


46, 


,9 


-03 


09 


24 


8 


.2 


6.7 


05/23 


18437-0216 


18 


46 


22, 


,7 


-02 


13 


24 


7 


.3 




05/23 


18440-0148 


18 


46 


36, 


.3 


-01 


45 


23 


8 


.3 




05/23 


18445-0222 


18 


47 


10, 


,8 


-02 


19 


06 


9 


.4 


5.3 


06/13 


18447-0229 


18 


47 


23, 


7 


-02 


25 


55 


8 


.2 


6.6 


05/30 


18449-0158 


18 


47 


35, 


,6 


-01 


55 


26 


8 


.7 


5.9 


06/13 


18454-0136 


18 


48 


03, 


,7 


-01 


33 


23 


11 


.9 


2.7 


05/30 


18454-0158 


18 


48 


01, 


,3 


-01 


54 


49 


5 


.6 




06/13 


18460-0307 


18 


48 


39, 


,2 


-03 


03 


53 


9 


.5 


5.2 


05/30 


18470-0044 


18 


49 


36, 


,7 


-00 


41 


05 


8 


.2 




06/13 


18472-0022 


18 


49 


50, 


,7 


-00 


19 


09 


11 


.1 


3.2 


05/30 


18488+0000 


18 


51 


24, 


.8 


+00 


04 


18 


8 


.9 


5.4 


06/13 


18521+0134 


18 


54 


40, 


,8 


+01 


38 


01 


9 


.0 


5 


05/30 


18530+0215 


18 


55 


34, 


,2 


+02 


19 


08 


8 


.7 


5.1 


06/13 


18540+0220 


18 


56 


35, 


,6 


+02 


24 


54 


10 


.6 


3.3 


05/30 


18553+0414 


18 


57 


53, 


,0 


+04 


18 


06 


12 


.9 


0.6 


06/19 


18566+0408 


18 


59 


09, 


,9 


+04 


12 


14 


6 


.7 




05/30 


19012+0536 


19 


03 


45, 


,1 


+05 


40 


40 


8 


.6 


4.6 


05/23 


19035+0641 


19 


06 


01, 


.1 


+06 


46 


35 


2 


.2 




05/23 


19074+0752 


19 


09 


53, 


.3 


+07 


57 


22 


8 


.9 


3.7 


05/23 


19175+1357 


19 


19 


49, 


.1 


+14 


02 


46 


10 


.6 




05/23 


19217+1651 


19 


23 


58, 


,8 


+ 16 


57 


36 


10 


.5 




05/23 


19220+1432 


19 


24 


19, 


,7 


+ 14 


38 


03 


5 


.5 




06/13 


19266+1745 


19 


28 


54, 


,0 


+ 17 


51 


56 


10 


.0 


0.3 


05/30 


19282+1814 


19 


30 


28, 


,1 


+ 18 


20 


53 


8 


.2 


1.9 


05/30 


19403+2258 


19 


42 


27, 


,2 


+23 


05 


12 


6 


.3 


2.4 


06/13 


19410+2336 


19 


43 


11, 


,6 


+23 


44 


06 


6 


.4 


2.1 


05/30 


19411+2306 


19 


43 


18, 


.1 


+23 


13 


59 


5, 


.8 


2.9 


06/13 


19413+2332 


19 


43 


29, 


.0 


+23 


40 


04 


6 


.8 


1.8 


05/30 


19471+2641 


19 


49 


09, 


.9 


+26 


48 


51 


2 


.4 




06/13 


20051+3435 


20 


07 


03, 


,8 


+34 


44 


35 


3 


.7 


1.6 


05/11 


20081+2720 


20 


10 


11, 


,5 


+27 


29 


06 





.7 




05/11, 


20126+4104 


20 


14 


26, 


,0 


+41 


13 


31 


1 


.7 




05/11 


20205+3948 


20 


22 


22, 


,0 


+39 


58 


05 


4 


.5 




05/11, 


20216+4107 


20 


23 


23, 


,8 


+41 


17 


40 


1 


.7 




05/11 


20293+3952 


20 


31 


10, 


,7 


+40 


03 


10 


2 


.0 


1.3 


05/11 


20319+3958 


20 


33 


49, 


.4 


+40 


08 


45 


1 


.6 




05/11 


20332+4124 


20 


35 


00, 


.5 


+41 


34 


48 


3 


.9 




05/23 


20343+4129 


20 


36 


07, 


.1 


+41 


40 


01 


1 


.4 




05/23 


22134+5834 


22 


15 


09, 


.1 


+58 


49 


09 


2 


.6 




05/06 


22551+6221 


22 


57 


05, 


,2 


+62 


37 


44 





.7 




05/06, 


22570+5912 


22 


59 


06, 


,5 


+59 


28 


28 


5 


.1 




05/06, 


23033+5951 


23 


05 


25, 


,2 


+60 


08 


11 


3 


.5 




05/06 


23139+5939 


23 


16 


09, 


,3 


+59 


55 


23 


4 


.8 




05/06 


23151+5912 


23 


17 


21, 


,1 


+59 


28 


49 


5 


.7 




05/06 


23545+6508 


23 


57 


05, 


.2 


+65 


25 


11 


0, 


.8 




05/06, 



Table 1. Positions of the IRAS point sources satisifying 
the Sridharan ct al. (2002) criteria, processed to J2000 co- 
ordinates, alongside the kinematic distance of the IRAS source 
and date(s) of observation. All sources were observed during 
the spring and summer of 2000. All distances are taken from 
Sridharan et al. 2002 with the exception of IRAS 18449-0158, 
for which we derive the kinematic distance using Vlsr from 
the CS(2-1) observations of Bronfman et al.(1996). Candidates 
with only the far kinematic distance listed have had their dis- 
tance uncertainty resolved. 



^evel of 0.03 Jy beam~^ and 0.69 Jy beam"^ was found 
_at 850 jum and 450 /zm respectively. The RMS noise level 
measured in each jigglcmap is listed in Table 2 as the un- 
certainty in the peak flux. 

Clumps were identified using object detection routines 
in the software package GAIA (Chipperfield & Draper 
2001). We define a positive detection as a group of pix- 
els subtending at least the area of the JCMT beam with 
emission above a 3-a level, where a is the RMS noise level 
of the jigglcmap. The vaUdity of each detection was also 
confirmed manually. Clumps not quite bright enough to 
be automatically detected were examined, and if deemed 
worthy of inclusion, added to the list of detections. These 
05/30 lower sigma detections are labelled by a note in Table 2. 

We list the peak fiux per beam and the integrated fiux 
for each detection. The peak flux per beam gives the peak 
flux level averaged in a 14.4" beam for 850 fim maps and 
in an 8.0" beam for 450 /xm maps, while the integrated flux 
of a detection measures the total flux inside an isophote 
tracing the S-cr RMS noise level around the detection. We 
quote the position of each detection as the location of peak 
emission, not as the centroid of the 3-a isophote. 

Calibrating the integrated flux of a detection required 
additional consideration, as the JCMT beam structure is 
complex (Figure 1), so the number of detector counts re- 
covered within an aperture is also a function of aperture 
size;. We quantified the extent of this relationship using 
maps of Uranus (which we consider a point source), cali- 
brating detector counts recovered inside circular apertures 
of increasing radius. We did not include the small number 
of non-planetary fiux calibrators in the calibration of in- 
tegrated fiux. The resulting function measures increasing 
counts with aperture size, asymptotically reaching max- 
imum counts once the aperture has expanded to encom- 
pass the JCMT beam and its primary error beam. For 
each detection, we then converted n counts recovered in- 
05/23 isophote of area A to Janskys by multiplying n by 

the counts-to-Jy conversion factor derived from a circular 
aperture of equivalent area. 

Our observations were performed on seven nights over 
a period of three months. Despite the protracted nature 
of our observations, a comparison of the counts-to-Jy con- 
05/23 version factor calculated for each night showed it usually 
05/23 ].gjjja,jjjg(j consistent with the published JCMT response.^ 
Where the conversion factor appeared inconsistent and 
^'^ other recent fiux calibrator maps were available, we 

assumed a conversion factor equal to the mean value for 

our run. A comparison with the independent 1.2 mm con- 
tinuum observations of Beuther et al. (2002a) shows the 
data to be consistently calibrated (§6.1), and we estimate 
the absolute flux uncertainty to be ± 10% at 850 /xm and 
± 30% at 450 iJ,m. 



http://www.jach.hawaii.edu/JACpublic/JCMT/. . . 
. Continuum_observing/SCUBA/astronomy/calibration/gains.html 
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Fig. 1. Jiggle-maps of Uranus measured at 850 /xm (left-hand 
image) and at 450 /xm (right-hand image). Contours were cho- 
sen to highlight the JCMT beam structure, and are drawn 
every 0.8 magnitudes below 65 Jy for the 850 jj,m image and 
every 0.8 magnitudes below 150 Jy for the 450 fim map. The 
primary error beam can clearly be seen as a ring encircling the 
main beam, containing ~ 9% and ~ 25% of the total flux at 
850 /Ltm and 450 /nm respectively. 



2.1. Telescope and reduction artefacts 

We masked a small number of consistently noisy 850 /um 
bolometers during data reduction; these bad bolometers 
usually fell in a region of background sky or faint ex- 
tended emission. Unfortunately, the secondary source seen 
towards IRAS 05358-1-3543 fell on a noisy 850 /xm bolome- 
ter, but wc consider the companion source an important 
feature of the jigglemap and so leave the bolometer un- 
masked. In addition, the map of IRAS 18553-1-0414 forms 
an isolated case that suffers from an imusually large num- 
ber of bad 850 /xm bolometers; we still include this data 
as the 450 fim map reveals the majority of flux has been 
recovered by good bolometers. 

Our observations used a 120" chop to sky to mea- 
sure and remove the background emission. However, in 
crowded regions, the 120" chop-throw sometimes points 
the telescope towards an occupied region of sky rather 
than an empty field. When this occurs, emission from ob- 
jects in the sky reference beam is subtracted from the tar- 
get field emission, resulting in negative images of clumps 
seen towards the reference position superimposed onto the 
final map. Some of our maps c;ontain these artefacts, which 
are usually seen away from regions of interest (eg. IRAS 
18151-1208, IRAS 18431-0312), but chopping onto emis- 
sion altered the map of IRAS 18454-0158 to such an extent 
that no reliable measurement was possible, and this source 
was removed from our analysis. 

We occasionally observed additional jigglemaps offset 
from the target position to map fields with emission con- 
tinuing outside the ~ 120" SCUBA field of view. These 
additional maps were calibrated as individual jigglemaps 
before they were combined into a mosaic, weighting the 
contribution of each map to intersecting areas by 1/cr^, 
where a is the RMS noise level in the map. Detections 
within the mosaic are still defined as clumps with emis- 
sion above a 3-a limit over an area the size of the JCMT 



beam, but using the RMS local to the section of mosaic 
being measured. 

3. Results 

Submillimetre emission at 850 /im was detected towards 
all the IRAS sources in our sample, although not all 
sources were bright enough to be detected above the in- 
creased background emission at 450 /im. When sources 
were bright enough to be detected at 450 /xm, the increased 
resolution of the 450 /xm observations sometimes resolved 
additional peaks within the area of a single 850 /im detec- 
tion. A presentation of the reduced jigglemaps, calibrated 
in Jy beam~^, alongside maps of a, the spectral index of 
the dust emission (detailed in §4.1), can be found in Figure 
2^. The position and flux measured for each detection is 
hsted in Table 2. 

The target sources IRAS 19266-1-1745 and IRAS 
18553-1-0414 displayed a gas-|-dust mass incompatible with 
the luminosity of the driving protostar, unless these 
sources are projected to the far kinematic distance (§6.2). 
Therefore, we reject the near kinematic distance for these 
objects and consider them resolved to the far kinematic 
distance for all subsequent analysis. 

3.1. Companion Clump Fraction 

The majority of candidate HMPOs were detected as com- 
panionless, compact and approximately spherically sym- 
metric submm clumps (eg. IRAS 05553-1-1631), although a 
signiflcant number exhibited submm nebulosity (eg. IRAS 
18566-1-0408), appeared in fllaments of emission (eg. IRAS 
18437-0216) and/or existed with multiple additional de- 
tections within the field of view (eg. IRAS 23545-^6508). 
While our sample was explicitly constructed to consist of 
isolated, companionless candidates away from sources of 
confusion, we found only 38 of the 68 target IRAS fields 
contained a single, companionless clump. The remaining 
IRAS fields contained more than one submm clump, usu- 
ally two detections, with the mosaic map towards IRAS 
18089-1732 containing the most companions, where 5 sep- 
arate clumps were resolved. This demonstrates the diffi- 
culty in locating truly isolated candidate HMPOs. Wc can 
characterise the multiplicity of our detections by calculat- 
ing the companion clump fraction (CCF), expressed by 
the formula 



CCF ■ 



B + 2T + 3Q + 4:P 
S + B + T + Q + P' 



(1) 



where S, B, T, Q and P are the number of single, binary, 
triple, quadruple and quintuplet clumps in our sample. If 
all clumps in our sample were solitary clumps, the CCF 
would be 0.0, while if all clumps had one companion the 
CCF would be 1.0. 



^ To meet size constraints. Figure 2 has been trun- 
cated to a single page; the complete figure displaying 
all 68 maps can be found in the preprint available at 
http://saturn.phy. umist.ac.uk:8000/~tjm/hmpol.ps.gz 
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Fig. 2. Maps of Saso, S450, and the a distribution in the left-hand, centre, and right-hand panels respectively. Greyscale limits 
are chosen to emphasise flux levels between -2 a and +7 a, where cr is the RMS noise level in the map. Contours trace the 
intensity in units of Jy beam~^, using the step-size listed in the bottom left-hand corner of each map. The first contour is drawn 
at the first step above zero Jy beam^^ unless the index is marked with an astcrix: this signifies that an additional contour is 
plotted at half a step above zero Jy beam~^. Triangular symbols plot the location of MSX point sources in each field of view. 
The spectral index of the dust emission (a) is plotted in the right-hand greyscale maps: a is masked outside the boundary of the 
first 850 /xm and 450 /xm contours. Contours on the a map directly mirror those drawn on the 850 /xm submillimetre emission 
map. 



Table 2. The position and measured flux of the detections resolved by our JCMT observations. Positions are measured from the 450 |a.m jigglemaps wherever 
possible. Unless otherwise stated, 450 \im detections without corresponding 850 \im detections arise from the increased resolution of the 450 |J.m observations. 
^850 |xm emission falls on a noisy bolometer. 

^2-0- detection at 450 (xm, where a is the RMS noise level in the map. 
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The CCF for our sample is 0.65 ±0.1, wliere the uncer- 
tainty comes from ^/N counting statistics. In reality, the 
absolute value of the CCF and the quoted uncertainty are 
both lower limits, as they are calculated assuming we have 
detected, and are uniformly sensitive, to all companions. 
This is not the case, as our limited angular resolution pre- 
cludes the detection of companions closer than around a 
beamsize, plus the finite field of view means companions of 
greater than ^ 60" separation (assuming a clump central 
in the jigglemap) will not be detected. Additionally, when 
coupled with the large difference in projected distance (the 
most distant sources being more than 15 times further 
away than the closest sources), our companion mass sen- 
sitivity also bears a dependence on distance. 

Wc examined the effect of different distance projec- 
tions by sorting our candidate HMPOs into four bins, 
containing sources < 2 kpc, 2 to 4 kpc, 4 to 8 kpc and 
> 8 kpc distant, respectively. The CCF of these subsam- 
ples remains remarkably consistent, each section in agree- 
ment with the full sample CCF within the uncertainty lim- 
its. This is true regardless of whether distance-unresolved 
sources are projected to the near or far kinematic distance, 
with the exception of the < 2 kpc bin projected to the 
far kinematic distance, and suggests c;lunips have a simi- 
lar number of companions over a wide range of distance 
scales. 

The expression of multiplicity given in Equation 1 is 
usually used as a diagnostic of more evolved stars, in par- 
ticular to quantify the number of companions a low-mass 
star is born with (eg. Beck et al. 2003, Patience et al. 
2002), whereas in this study the CCF can be interpreted 
as the likelihood of finding additional potentially star- 
forming clumps when observing Galactic HMPOs iden- 
tified by a similar flux-limited criteria. While it remains 
difficult to constrain the statistics of such a disparate sam- 
ple, the CCF does emphasise that most clumps do not 
form in isolation, and that a single IRAS detection is usu- 
ally resolved into several submm clumps. New samples of 
protostars comparable to our sample are hard to compile, 
but the strong likelihood of detecting additional clumps 
in the locality of our sample suggests that wide-field sur- 
veys towards existing high-mass protostars may also be a 
productive way of locating new protostellar candidates. 

3.2. Clump morphology statistics 

Forming an unbiased statistical analysis of the morpho- 
logical features is difficult, as distinguishing features are 
mainly found in the appearance of low-level extended 
emission. As such, measuring the FWHM of detections 
is of limited use, as it is not sensitive to the faint emis- 
sion features we wish to characterise. Instead, we formed a 
simple statistic that indicates how much mass lies outside 
the central beam by measuring Y, the ratio of integrated 
flux (Sint) to peak flux (Speak) at 850 fim. We measure 
at 850 /Ltm because of higher signal-to-noise than in the 
corresponding 450 i^m maps. For a point source, Y equals 
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5 10 15 

Sipt/Speok (14.4" beams) 

Fig. 3. A histogram of the Sint / Speak flux ratio for 'confirmed' 
isolated sources (sample A: shaded bars) and sources poten- 
tially with companions (sample B: outline). 



unity. The value of Y for each detection is listed in column 
9 of Table 2. 

In Figure 3 we plot the Y distribution for our submm 

detections, dividing our sample into two groups: subsam- 
ple A, containing detections with a high confidence of 
being solitary cores (having just one detection within an 
IRAS, MSX, SCUBA and IRAM field of view) and sub- 
sample B, the remainder. While the solitary detections 
in sample A do not display the extended distribution tail 
seen in subsample B, we see both groups peak at an in- 
tensity ratio of ~ 3, which in light of the factor 15 range 
in distance suggests that the envelope structures may be 
scale- free. 

The large scatter in Figure 4 means we do not find any 
strong correlation with distance at the near distance pro- 
jection (which is the most likely projection for our sam- 
ple), although considering the distance-resolved sources 
alone does reveal a trend beyond d = 4 kpc, where 
Y (X (i""-^. This fall-off does not reveal a physical change, 
but reflects the diminishing level of integrated emission 
as the 3-ij isophote encloses less of the envelope for more 
distant sources. 

We must qualify a number of uncertainties that could 
affect the distribution of Y, not least our variable sensitiv- 
ity to additional embedded sources. While the large-scale 
envelope structures appear to be scale-free, on the small 
scale there are indications that we are still undersampling 
the number of companions separated by less than a beam 
width. A number of apparently single detections at 850 
yum are barely resolved as multiple sources at 450 /xm (eg. 
IRAS 05490-1-2658); a reminder that further clustering on 
size scales less than a 450 fim beam width may also be 
present. Overall, it is inevitable that with limited resolu- 
tion we misclassify some multiple cores as solitary detec- 
tions, a point demonstrated in Figure 5, where all sources 
(bar one) with 850 /zm flux ratios larger than 7 are resolved 
as muhiple detections at 450 fim. IRAS 22551-^6221 pro- 
vides the most visible demonstration of this effect, where 
the high 850 /im flux ratio arises from the inclusion of 
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Fig. 4. A plot of y = Sint/ Speak, tho ratio of 850 /^m integrated flux to peak 850 jira flux measured iu a 14.4" beam, against the 
kinematic distance of each detection. Distance resolved sources are plotted by filled circles, while distance unresolved sources are 
projected to near and far kinematic distances and plotted with open rectangles and triangles. The curve in the distance-resolved 
plot displays a power law of the form Si„t/ Speak oc d~°'*. 




peak 



Fig. 5. A comparison of the peak flux and integrated flux 

of each 850 /im detection. The solid line traces the Y = 
Sint/ Speak = 3 distribution peak found in Figure 3, while the 
dotted line traces the ratio Y = 7. 



flux from a bright neighbouring source that is only fully 
resolved at 450 /um. 

We also tend to overestimate the flux of multiple de- 
tections, as the elliptical apertures used for photometry 
could include emission shared with a companion source. 
Although the intersection of apertures around adjacent 
components was minimised where possible, it remains a 
potential cause of uncertainty. Finally, the ratio for ex- 
tended sources is likely to be a lower limit, as emission 
from a large, extended envelopes is more likely to project 
emission onto a noisy bolometer, and flux incident on these 
noisy bolometers is masked during jigglemap reduction. As 
a result, the quoted integrated emission is a lower limit, 
and the flux ratio is underestimated. 

To conclude, while these concerns affect the quanta- 
tive results, qualitatively we still observe that a significant 
fraction of the total mass lies outside the central 'core' at 
this stage of evolution. 

3.3. Surface Density 

Stars generally form in association with other stars, and 
the spatial distribution of these groups of stars can provide 




■"0.1 



10 



separation (pc) 



Fig. 6. MSDC for single power-law distribution models (light 
dotted lines), with p(N< r) oc from 7=-2 (uppermost dotted 

line) to 7=3.5 (bottommost dotted line) with Sj = 0.5. The 
square symbols represent the observed MCSD (multiplied by 
100), with a thick dashed line plotting the line of best flt for the 
observed MSDC above the break-point. The best fit power-law 
distribution (7=-0.75) is plotted by a thick black line. 

information on how the natal molecular cloud fragmented. 
One way to probe the distribution of sources is to use 
the mean surface density of companion sources (MSDC), 
a method which has been successfully used to probe the 
transition from the formation of binary stars to star clus- 
ters in low mass star forming regions (eg. Gomez et al. 
1993; Larson 1995; Nakajima et al. 1998). 

The MSDC of our detections is shown in Figure 6. The 
MSDC was calculated by measuring the linear separation 
r of each detection to its companions. The separation of 
each companion pair was binned into annuli of separation 
r to r -I- Sr. The number of pairs A'" within each annulus 
was then divided by the area of the annulus and the total 
number of sources to give the MSDC, S(r), as 



S(r) = N/{2T:rdrN^) 



(2) 



Above 0.4 pc, the MSDC of our detections has a mea- 
sured gradient 7 of -1.7, which is roughly halfway between 
the power-law indices of binary pre-main-sequence popu- 
lations (where 7 —0.5; Nakajima et al. 1998) and that 
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of more distant companions (where 7 ~ —2.2; Nakajima 
et al. 1998). The MSDC we observe appears to turn over 
below 0.4 pc, but the validity of this turnover is ques- 
tionable as we have few measurements in this region. The 
clump MSDC break is also quite distinct from the stellar 
MSDC power-law break found at ^ 0.04 pc (Gomez ct 
al. 1993; Larson 1995, Nakajima et al. 1998), below which 
the MSDC steepens from the inclusion of close binaries, 
as the stellar MSDC breakpoint occurs on much smaller 
scales than are detectable by our survey. 

To identify the space density distribution consistent 
with the clump MSDC, wc modelled a munbcr of systems, 
each containing a collection of 10'^ sources randomly dis- 
tributed according to a number of power-law space den- 
sity distributions. The projected MSDC of these simula- 
tions are also shown in Figure 6. For separations above 0.4 
pc, the MSDC power-law slope of -1.7 most closely cor- 
responds to a space density distribution with the number 
density per unit volume Nyoi(r) oc 7-"°-^^. 

We must consider the MSDC statistic cautiously, for 
the MSDC is constructed from observations of widely dif- 
ferent companion seperation sensitivities and companion 
flux sensitivities, in a similar vein to the CCF (§3.1). 
Complications arise from the wide range of projected dis- 
tances to our HMPO candidates: specifically, our observa- 
tions are sensitive to angular separations from around a 
beamwidth up to the upper limit of a 120" field of view, 
but as most HMPO candidates in our sample are less than 
4 kpc away, there are only a very small range of uniformly 
sampled linear separations for our candidates. Towards 
more distant IRAS fields, we are uniformly sensitive to 
a larger range of linear separations, but we then suffer 
from fewer measurements and from a reduced sensitivity 
to close companions. 

We examined the significance of variable sensitivity us- 
ing a procedure similar to that used for the CCF (§3.1), 
comparing the MSDC slope of groups of sources with sim- 
ilar distances, finding the slope of each MSDC segment 
agrees with the MSDC of the whole sample within the 
uncertainty limits. We suggest the MSDC as calculated 
provides at least a basic estimate of the clustering proper- 
ties of these clumps. Ultimately, the scarcity of high-mass 
protostars means there will always be a large range of dis- 
tances in samples of HMPOs, and we may never be able 
to construct a set of uniformly sampled observations to 
the extent possible with low-mass protostars. 

3.4. Clump positions 

For each candidate HMPO, the telescope was pointed so 
that the IRAS point source was central in the field of 
view, hence the location of the IRAS source in each map 
is given implicitly by the map centre. We indicate the posi- 
tion of neighbouring MSX sources in the field of view with 
triangular symbols in Figure 2, finding that many IRAS 
sources and MSX detections are roughly coincident with 
the submm clumps, but in agreement with SBSMW we 



also find that some submm detections and IR detections 
cannot be coincident within the positional uncertainties of 
our survey and the IRAS/MSX surveys. For example, to- 
wards IRAS 23545-1-6508 there are two MSX point sources 
found within 30" of the two submm clumps, an ofii'set 
greater than the expected absolute uncertainty, plus the 
MSX point sources also have smaller relative separation 
between components than the corresponding SCUBA de- 
tections. Sources with large SCUBA/MSX offsets do not 
appear to have /3 or Ss5o characteristics different to more 
coincident detections. 

Considering each clump as a protocluster may explain 
displaced IRAS/MSX and submm detections, as addi- 
tional stars embedded in the less dense, more transparent 
outer reaches of the envelope will not encounter the same 
degree of opacity, providing a mechanism for shorter wave- 
length photons to pass. To examine this possiblity and 
resolve whether MSX and SCUBA detections trace the 
same body of material will require further high-resolution 
IR observations. 

4. Dust optical depth 

The dust optical depth r can be determined from the flux 
density using the expression 



where B^(T) is the Planck function at temperature T, and 
il is the solid angle subtended by the telescope beam, all 
measured at frequency ly. This relationship assumes the 
emission fills the telescope beam, which may not be the 
case, so the derived value of t^, must be considered the 
beam-averaged optical depth. 

To calculate the dust optical depth for our detections 
using Equation 3 we assumed dust temperatures equal to 
the SBSMW cold-component dust temperatures. In the 
SBSMW study, the spectral energy distribution (SED) of 
each IRAS source was successfully modelled as a compos- 
ite of two greybodies: one greybody representing a cold 
dust component, accounting for the > 60 /xm flux, while 
a separate hot dust component contributes the majority 
of near-IR flux. As SCUBA is only sensitive to emission 
from the cold dust greybody, we set T eqiial to the tem- 
perature of the cold component {Ted) as given in Table 1 
of SBSMW. 

The beam-averaged optical depth and other param- 
eters derived from the flux density are listed in Table 3. 
While our sources consist of very dense clumps, the beam- 
averaged 450 ^m and 850 /xm optical depths show that 
they are usually optically thin at submm wavelengths. The 
optical depth at 850 ^m spans almost two orders of mag- 
nitude, from 10"'^ < T850 < 10~^. At 450 ^m, T450 is found 
within the range 2 x 10~^ < T450 < 10~^ with three excep- 
tions: detection #8 (IRAS 18089-1732) appears optically 
thick at 450 //m with r45o = 1.1 but this is by far the 
brightest detection of our survey. Two other detections 
(#16: IRAS 18182-1433 and #19: IRAS 18264-1152) have 



Table 3. Derived parameters of the submni detections resolved by this survey. The mass of each clump is calculated from the 850 /um integrated flux value, 
assuming the dust grains have thin ice mantles, and using a 100:1 gas-to-dust ratio. The column density refers to the total gas column density (ie. n(H-|-H2)). 
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high 5*450, leading to higher T450 than the majority of the 
detections, but they remain with T450 < 1- 

Wc have modelled the submm emission seen towards 
our sample using a one-dimensional radiative-transfer 
code (details and results can be found in the companion 
to this paper: Williams et al. 2004). These models assume 
a fixed dust grain chemical composition, with a silicon-to- 
graphite ratio half that of the interstellar medium (Mathis 
et al. 1977) and a standard Draine & Lee (1984) dust grain 
size distribution. We used these model dust grains to pre- 
dict the optical depth of the cores as a function of wave- 
length, which when scaled to match the observed 850 /im 
and 450 /um optical depths gives an estimate of the wave- 
length at which the submm detections become optically 
thick (listed in column 7 of Table 3). We find that for the 
average clump, r > 1 for wavelengths shorter than ~90 
fim. 

4.1. Spectral index of the dust emission 

As noted in §4. the dust emission at submm wavelengths 
can be well represented by a greybody, with intensity vary- 
ing smoothly as a function of frequency. Having measured 
the intensity of dust emission at two submm frequencies, 
we can characterise the SED using a, the spectral index 
of the dust emission. This is defined 



HS1/S2) 



(4) 



ln(i/i/i/2) 

where Si and ^2 are flux densities at frequencies i>i and 
1^2 respectively. By dividing calibrated images measured 
at frequencies vi and 1^2, we can determine the spatial 
distribution of a and examine its relationship with the 
intensity of emission and density of gas and dust. 

However, first we must take into account the differ- 
ent JCMT response and beam patterns at 850 and 
450 /Ltm, for as seen in Figure 1, more flux lies outside 
the main beam at 450 fim compared to the response at 
850 /im. We accoimtcd for these differences by following 
the procedures defined in Hogerheijde & Sandcll (2000), 
normalizing the images to a common response before fi- 
nally determining a as the ratio of images. In detail, we 
described the JCMT beam at each wavelength as a su- 
perposition of three Gaussians, the parameters of which 
were found by a fit to the azimuthal average response to 
Uranus. The amplitude and FWHM of these components 
are listed in Table 4. We then deconvolved the 850/im 
and 450/Lim SCUBA images with the corresponding beam 
model, smoothing the deconvolved images with a single 
Gaussian to achieve a final, uniform, spatial resolution of 
15.0" before forming a as given in Equation 4. 

The spatial distribution of a can be seen to the right 
of the submm emission maps in Figure 2; a is blanked for 
emission lying outside the first 850 /im and 450 /xm con- 
tours, so only the intersecting area of the maps with good 
signal to noise is both displayed and analysed. In Table 3 
we list the spectral index both in terms of the index at the 
position of peak 850 /xm flux {apeak', column 13) and the 



850 




450 




Relative 


FWHM 


Relative 


FWHM 




I, ) 




I, ) 


0.93 


14.4 


0.81 


8.0 


0.05 


44.3 


0.17 


24.7 


0.02 


62.9 


0.02 


71.3 



Table 4. Parameters used for a three Gaussian component 
description of the JCMT beam. 



mean value of the a map {amean, column 12). Unless we 
state otherwise, the spectral index is discussed in terms of 
apeak, the index at the location of peak submm flux. Also, 
we do not include detection #13 (IRAS 18102-1800) in the 
analysis, as it appears an outlier with significantly lower 
a than any other detection; we suspect this is caused by 
the suspiciously weak 450 /im emission seen towards this 
source, most likely due to a transient telescope or calibra- 
tion error as the 450 /xm flux appears inconsistent with 
the 100 /zm, 850 /im and 1 mm flux constraints. 

The spectral index measured at the position of 
peak 850 fxm emission varies from apeak = 11 (IRAS 
19411-F2306) to apeak = 3.7 (IRAS 20293-^3952), though 
overall this index is fairly uniform with a sample mean of 
apeak = 2.6±0.4. Averaging all spectral index data around 
a detection slightly reduces the statistical variability, so 
that amean ranges from 1.1 to 3.5, and the sample mean 
falls to amean = 2.4 ± 0.3. As imphcd by these statistics, 
the majority of detections display an a distribution that 
peaks towards the location of maximum 850 fj.m emission, 
although some sources display an anticorrelation with in- 
tensity. 

Cross-sections of the a distributions towards IRAS 
05358-^3543 and IRAS 05490-^2658 are presented in 
Figure 7. Although these sources display very different 
a morphologies, we note that the positive and negative 
features seen towards the location of peak submm emis- 
sion are roughly comparable in depth and width. As these 
sources are at roughly the same kinematic distance, the 
cause of these features could potentially occur on a simi- 
lar spatial scale. 

4.1.1. Peaked a distributions 

The spectral index of the dust emission depends on a com- 
bination of the beam-averaged values of dust temperature, 
opacity, and spectral index of the dust opacity (/3, deflned 
in the sense t oc v^). As a result, there are three mecha- 
nisms which may explain the spatial distributions we ob- 
served: 

1. Temperature variations through the dust envelope. 

2. Emission originating from optically thick regions. 

3. Changes in the composition of the dust grains them- 
selves. 
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Fig. 7. A plot of the typical cross-sections seen for peaked a morphologies (IRAS 05358+3543; left-hand plots) and negative- 
dip morphologies (IRAS 05490+2658; right-hand plot). The upper row displays the orientation of the cross-sections, while the 
bottom row displays the a index measured along the labelled cut, using a pixel scale of 3" per pixel. 



The only currently available estimate of dust temper- 
ature was derived by SBSMW using greybody fits to the 
SED at long wavelengths. The majority of dust cores 
have a temperature similar to the mean of the sample 
(7d„st =44 K), although a number of cores are associ- 
ated with higher temperatures: IRAS 18440-0148 with 
Tdust=97 K, IRAS 20319+3958 with Tdust=73 K, and 
IRAS 23151+5912 with Tdust=GS K are the most promi- 
nent higher temperature cores. However, with only a single 
dust temperature estimate for each source, the magnitude 
of any dust temperature gradient across the protostellar 
envelopes remains unknown, which leaves the contribution 
of any temperature-dependent mechanism to variations in 
a unclear. 

Single-dish NII3 observations have also been con- 
ducted towards our sample, tracing gas within the cooler, 
extended envelope (SBSMW). From these observations, 
SBSMW found a mean temperature of Tnh3—^Q K, 
around 25 K lower than the dust temperature in an av- 
erage core. However, a core containing warm dust and 
an extended envelope characterised by cool gas does not 



prove the existence of a temperature gradient, for it is 
very difficult to make the gas temperature close to that of 
warm dust, even with the high densities {N^ ~ 10^ cm~^) 
seen towards typical protostellar candidates (Goldsmith 
et al. 1997). However, CH3OH and CH3CN molecular 
tracers, pointing to high temperature, high density re- 
gions and indicating the presence of a hot core, have 
been detected towards a number of our candidate HMPOs 
(SBSMW). These detections imply that there are regions 
within the clumps of much higher temperature than the 
beam- averaged dust temperature alone would suggest, so 
we expect a strong temperature gradient must be present 
towards at least some of our sources. 

We examined the significance of a temperature gradi- 
ent by forming simulated a maps using the 450 /im and 
850 /im continuum images created by our best-fit radia- 
tive transfer models (Williams et al. 2004). These models 
assume a single luminous protostar embedded in a dense, 
dusty envelope, and form excellent fits to the observed 
emission while maintaining constant dust grain character- 
istics (ie. opacity and variable grain composition are not 
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a factor in the simulated a map). In general, our best 
fit models suggest the presence of dust envelopes with 
temperatures around 300 — 500 K at the inner bound- 
ary, falling to around 10 — 15 K at the outer boundary. 
With a temperature gradient as the only factor, the sim- 
ulated a distribution peaks towards the hottest, densest, 
most central regions. These centrally peaked a morpholo- 
gies are similar to those seen towards the majority of our 
sample (eg. IRAS 18247-1147; IRAS 18306-0835), suggest- 
ing the observed a features are dominated by temperature 
gradients across the envelope. This result emphasises that 
we must know the spatial temperature distribution of the 
clumps if we are to refine our investigation and accurately 
quantify the contribution of other factors towards these 
sources. 

4.1.2. a-dip distributions 

On the other hand, peaked a distributions are not the 
only morphology observed: the IRAS sources 05490-1-2658, 
18290-0924, 18530-1-0215, 19413-^2332 and 20051-^3435 
form notable exceptions where a falls towards the loca- 
tion of maximum 850 /zm emission. In terms of tempera- 
ture gradients, these a morphologies run counter to the a 
distribution expected for an internally heated core: if no 
other factors are involved, they imply that the inner core 
must be cooler than the surrounding envelope. But is the 
formation of a hot envelope and cool inner core a realis- 
tic possibility? To form a typical negative dip morphology 
with an a valley depth of Aa = —0.3, while maintaining 
an average dust temperature of 44 K. would require an 
inner core temperature of around 26 K with a surround- 
ing envelope of 60 K. It would be hard to explain such 
low inner temperatures in the presence of large, luminous 
protostars when 'hot-cores' associated with typical pre- 
UCHII protostars have temperatures of > 100 K (Kurtz 
ct al. 2000). Alternatively, external heating could warm 
the outer layers of the clump relative to the inner core, 
but the interstellar radiation field alone is not capable of 
heating such dense dust to such high temperatures, and 
while nearby luminous stars could conceivably heat the ex- 
terior to higher temperatures our radiative transfer mod- 
elling shows the submm emission profiles arc well matched 
by low temperatures (~ 10 — 15 K) at the external enve- 
lope boundary (Williams et al. 2004). Furthermore, low 
core temperatures may prcchide the very formation of a 
massive protostar, as the Jeans mass becomes much lower 
within cooler cores, suggesting that a series of lower mass 
protostars would form instead. This does not exclude the 
possiblity that a massive star could form through the co- 
alescence of discrete low-mass protostars, but considering 
the weight of evidence it is hard to envisage how 'cool 
cores' may cause the observed variations. 

An alternative explanation is that these a-dip cores 
are optically thick. This possibility is unlikely, as §4 shows 
that all cores (with the exception of Source #8, the main 
component towards IRAS 18089-1732) are optically thin 



even at 450 /um. Even so, we recognise that this state- 
ment is based on the beam-averaged values, and there 
may be much denser, optically thick regions present on 
scales smaller than our observations can probe. For exam- 
ple, a circumstellar disk would lead to a large density con- 
centration in the very inner envelope, but would remain 
unresolved by our observations. Then again, the presence 
of circumstellar disks has been confirmed towards IRAS 
20126+4104 (Cesaroni et al. 1997) and IRAS 05553+1631 
(Shepherd & Kurtz 1999), and the a distribution towards 
these sources remains strongly peaked. Overall, we con- 
clude that optically thick regions do not significantly affect 
the a distribution at the spatial resolution of our measure- 
ments. 

Finally, variations in the properties of the dust-grains 

themselves could help explain the trends. The optical 
properties of dust grains can be quantified by /3, the spec- 
tral index of the dust opacity. This is often a quantity of 
interest as it may give information on the composition and 
evolutionary history of dust grains within the envelope. 
There are many models that predict (3 for different grain 
characteristics, and the majority of grain compositions re- 
sult in a spectral index of /3 ~ 1.5 — 2 (eg. Gezari, Joyce & 
Simon 1973; Draine & Lee 1984; Kriigel & Siebenmorgen 
1994), although it may range from /3 ^ 1 (Mathis & 
Whiffen 1989) up to /3 = 3 (Aannestad 1975). To observe 
an a-dip morphology, dust grains within the central core 
must be of lower P than grains in the surrounding enve- 
lope. Low P and /3 distributions that fall towards regions 
of high density are usually attributed to grain growth in 
these dense, innermost regions (eg. Mannings & Emerson 
1994; Beckwith & Sargent 1991; Goldsmith et al. 1997), 
and our observed a-dip morphologies are generally consis- 
tent with this grain growth interpretation, as a (and thus 
P) fall preferentially towards the centre of the dense cores 
we have observed. 

However, the models of Ossenkopf & Henning (1994) 
predict that /? will only change if the dust grains do not 
have ice mantles, which would require inner cores with 
dust temperatures > 100 K. The detection of CH3OII and 
CH3CN towards the a-dip detections IRAS 19413+2332 
and IRAS 18530+0215 signifies the presence of a hot core 
of sufficient temperature to melt ice mantles, thus permit- 
ting grain growth. However, the detection of CH3OH and 
CH3CN towards a large number of candidate protostars 
with centrally peaked a morphologies raises an interest- 
ing question: these molecular tracers imply a high central 
temperature - certainly high enough to melt ice mantles 
and permit changes in /3, so why are there no signs of 
grain growth? Why do the majority of these sources have 
positively peaked a morphologies? The strongly peaked 
a distributions we observe suggests that the temperature 
gradient effect outweighs any contribution from grain evo- 
lution. On the other hand, as it takes time for ice mantles 
to melt, perhaps these icy grains remain towards cores 
only recently heated, and perhaps these cores are younger 
than those associated with a-dip distributions. 
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Also, CH3OH and CH3CN have not been detected to- 
wards IRAS 05490+2658 and IRAS 18290-0924, suggest- 
ing a hot core has not formed, yet these objects are still 
found with a-dip morphologies. Clearly, neither temper- 
ature gradients or variable dust grain compoaition taken 
alone cannot fully explain the observed a morphologies. 
It is not clear that hot cores are strongly correlated with 
grain growth nor with a-dip morphologies, and accurate 
high-resolution measurements of the temperature of the 
clumps are vital if we are to determine the magnitude of 
grain growth towards our sample. 

4.2. P and signs of grain growth 

As we have an estimate of the temperature of dust grains 
within the cores, we can calculate the spectral index of 
the dust opacity using the equation 



/? (a + Aa) - 2, 



(5) 



where Aa is a Rayleigh-Jeans correction factor, neces- 
sary because the dust body temperatures are closer to 
the equivalent temperature Ti, of our observations than a 
Rayleigh-Jeans approximation would permit (where tem- 
perature at frequency v is = hv/k, giving Tg^o w 17 
K and T450 ~ 32 K). Assiiming the equivalent tempera- 
tures Ti^j and T,^^ for the observed frequencies and 1/2, 
the required correction can be expressed 



In 



Aa = 



cxp(T,2/Td„,t)-l 



oxp(T,i/Td„,t)-l 



ln(T^2/r^J 



- 1. 



(6) 



For example, with a median cold-component dust temper- 
ature of 43 K, the average correction is around -f 0.3. By 
taking our a maps and adding (Aa — 2) to each image, the 
a maps are transformed to display the spatial distribution 
of /3. However, as we only have an estimate of the dust tem- 
perature in the inner core (setting Tdusi = Td, the cold- 
component dust temperature from SBSMW), likewise the 
estimate of f3 obtained from these maps is only valid to- 
wards the innermost regions. 

The value of (3 at the location of peak 850 fxm submm 
emission is listed individually for each detection in Table 
3. The average grain opacity index for our sample is 
P = 0.9 ± 0.4. This index is smaller than for that seen 
towards high-mass stars associated with UCHII regions 
(Hunter 1997), and the /3 distribution is substantially 
shifted to lower indices compared to the more evolved ob- 
jects (Figure 8). Low /? is often associated with young, 
less evolved sources, which would point to further evidence 
that our sample of young high-mass stars are at an ear- 
lier stage of evolution than their UCHII counterparts. An 
inadequate Rayleigh-Jeans correction could raise our esti- 
mate of /3, but even if we have globally overestimated core 
temperatures by 20 K, /? would only to rise to ~ 1.2, still 
lower than /? towards more evolved UCHII sources. No re- 
lationship is found between a, (3 and distance (Figure 9), 
suggesting resolution is not an issue. 
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Fig. 8. A histogram of the Ppeak distributions for our sample 
of HMPOs (filled histogram) and for young high-mass stars 
associated UCHII regions (clear bars). The latter distribution 
is derived from the cold-component (3 values listed in Table 3.5 
in Hunter (1997). 



Figure 10 shows that despite the absence of any cor- 
relation between a and tsso or a, (3, tsso or Tdust and 
source distance, there is significant correlation between (3 
and r85o, albeit with a large scatter. Our estimate of (3 is 
dependent on the SBSMW dust temperature, which was 
calculated assuming (3=2. Had SBSMW used lower values 
of (3, they would have derived higher dust temperatures, 
which in turn both reduces the Rayleigh-Jeans correction 
and decreases the implied dust optical depth. As a result, 
I0W-/3 points in Figure 10 would move down (due to the 
smaller Raylciigh- Jeans correction) and to the left (due to 
decreased optical depth), a shift greatest for points to- 
wards the bottom-left quadrant as these objects would 
be associated with the largest dust temperature increase, 
hence strengthening the correlation between f3 and tssq. 

Low values of (3 are typically attributed to growth and 
evolution of the dust grains within dense, dusty regions, 
so the observed trend of low values of (3 at low optical 
depths is initially surprising. We would expect the dens- 
est, most massive clumps to undergo the most significant 
grain growth, resulting in lowest /? for these cores, whereas 
actually we seem to observe the opposite trend. In addi- 
tion, Ossenkopf & Henning (1994) found that dust grains 
within a protostellar core remain below the Rayleigh scat- 
tering size limit after 10^ years (the typical age expected 
for our sample; Behrend & Maeder 2001), with the dust 
grain opacity changing only by a factor of ~ 2 at 850 
/xm. However, the results in Figure 10 suggest a change of 
close to an order of magnitude in opacity: for values of (3 
around 0.5, rgso has a value of around 4 x 10^^, increasing 
to T850 ^ 3 X 10~^ for dust grains with /3 = 1.5. Only in 
the extremely dense and cold regions within circumstellar 
disks is it believed that a significant numbc;r of grains can 
grow beyond the Rayleigh limit, allowing the large shift 
in opacity (e.g. Schmitt et al. 1997). 
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Fig. 9. A plot of Qpeafe (left-hand plot) and Ppeak (centre and right-hand plot) against distance. Distance resolved sources are 
plotted by filled circles, while the square and triangular symbols represent distance unresolved sources projected to the near 
and far kinematic distance respectively. 
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Fig. 10. A plot of the relationship of Opeak (left-hand plot) and /3peofc (right-hand plot) to rsso, the optical depth at 850 /xm. 



These inconsistencies can be understood by consider- 
ing the inhomogeneous nature of our sample, and how 
the clumps we have observed vary in mass (§5) and most 
probably evolutionary status (SBSMW). Sources within 
the high T850, high /? quadrant of Figure 10 show consid- 
erable 'excess' mass compared to the mass of an equivalent 
luminosity main-sequence star (cf. Figure 11), suggesting 
these clumps could easily be forming protogroups or pro- 
toclusters (§5.1). In any case, the majority of dust grains 
within these high-mass clumps will not be intimately asso- 
ciated with the high- mass protostar, and will most likely 
remain outside the T > 100 K boundary necessary to melt 
ice mantles and allow large variations in (3. 

In contrast, the much lower mass of the low rgso, low 
/3 detections means these clumps may be forming solitary 
high-mass stars (solitary for high-mass stars being a rela- 
tive term, which we interpret as existing with only a small 
number of lower-mass stellar companions) , with fraction- 
ally much more dust lying inside the massive protostar's 
sphere of influence. This distinction becomes important 
when we consider that the T85o and P we observe are 
actually the optical depth and dust opacity index aver- 
aged along the line of sight. Within the largest clumps, 
low P grains in the vicinity of the high-mass protostar 
will be rendered less detectable, swamped by the higher 
/3 dust grains lying within the envelope of the proto- 



1 .5 



0.5 



□ 

□ □ 



□ • 



□ □ □ 

□ 



□ • 



10 



10" 



lass (I 



10^ 

^0) 



10 



Fig. 11. A comparison of detection mass and 3 at the location 
of peak 850 /nm emission {(3peak)- Sources with known distance 
are plotted by filled circles, while sources whose distance re- 
mains unresolved are projected to the near kinematic distance 
(square symbols) and the far kinematic distance (triangular 
symbols). 



group/cluster, whereas fractionally there will be many 
more low-/3 grains along the line of sight towards low mass, 
low T850 cores, making these evolved grains appear more 
prominently towards less massive cores. 
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Secondly, grain growth will occur in parallel with evo- 
lution of the high- mass protostar. The mass of material in 
the core will decrease as material either collapses to form 
stars or has been removed from the core by the action of 
stellar jets and winds. These mechanisms will act to re- 
duce the opacity towards more evolved cores by removing 
dust and gas. For example, the outflows of high-mass stars 
may disperse up to eight times the mass of material that 
falls onto the star (Churchwell 1997), allowing the dust 
opacity to change by a larger degree than that possible 
via grain growth alone. 

5. Mass 

The majority of photons we detect with SCUBA origi- 
nate from optically thin dust emission, so flux contours 
also trace the mass and column density within each field 
of view. The combined mass of the gas and dust can be 
calculated from the expression 

M = gS,£/>iMTdust), (7) 

where is the flux density, d is the distance to the 
source, is the absorption coefficient per unit mass of 
dust, g the gas-to-dust ratio and B,^{Tdust) represents the 
Planck function for a blackbody of temperature T^ust, 
all measured at frequency v. We adopt SBSMW kine- 
matic distances and cold-component dust temperatures, 
but as companion clumps do not have individual tem- 
perature measurements we have to assume all multiple 
detections in a field of view have the same dust tem- 
perature. We derive the absorption coefficient at 850 /im 
(k85o) from the evolutionary opacity models of Ossenkopf 
& Henning (1994), assuming an initial hydrogen number 
density of uh — 10® cm^'^, thin ice mantles and a forma- 
tion timescale of 10^ years to get Kgso = 1-54 x 10~^ cm^ 
g^^ and K450 = 5.23 x 10~^ cm^ g~^. For parity with other 
continuum observations, we assume a gas-to-dust ratio of 
100:1, and we list the derived clump masses in Table 3. 
However, note that the mass may be uncertain by a fac- 
tor of 2.5 or more, due to large variation in the gas-to-dust 
ratio (eg. Hildebrand 1983; McCutchcon ct al. 1995). 

Figure 12 shows the histogram of clump mass, where 
we see the majority of clumps have a mass of less than 
500 Mq, regardless of whether distance-unresolved cores 
are projected to the near kinematic distance or the far 
kinematic distance. The distance-resolved detections in 
our sample have a mean clump mass of around 350 
M0, although the median mass is less at around 100 
Mq. Assuming the near kinematic distance for distance- 
unresolved sources results in a mean clump mass of ~ 330 
Mq, with a median roughly half this value at 143 Mq, 
while projecting to the far kinematic distance results in a 
mean clump mass of 1120 Mq and a median mass of 460 

Mq. 

The mass of the distance-resolved and near distance- 
projected clumps are comparable to other continuum- 
derived mass estimates of similar massive protostars. For 
example, Mueller et al (2002) find an average clump mass 



of 209 Mq, while Molinari et al. (2002) measure M = 235 
Mq. These values support claims that natal clumps bear- 
ing massive stars may contain up to 100 times the mass 
of the most massive adult star that emerges (Churchwell 
1997). This is a very different scenario to that seen in the 
primarily low-mass star-forming region of Rho Ophiucus 
(Motte et al 1998), where a high fraction of the initial 
clump mass is seen to transfer onto the resulting low-mass 
protostars. 

Using the mass and column density relationships de- 
fined in Hildebrand (1983), we can also translate the mass 
of each detection to a beam-averaged gas column density. 
The values we derive are presented in Table 3, where we 
find an average H-I-H2 column density of 9 x 10^^ cm~^. 
With an average clump diameter of 0.6 pc and projected 
distance of 4 kpc, this translates to a mean hydrogen num- 
ber density of 3 x 10^ cm~^ through the clump. 

5.1. Cumulative mass spectrum 

Figure 13 presents the cumulative mass spectrum of our 

850 /im detections, distance-limited to IRAS fields less 
than 5 kpc distant to increase the region of complete sam- 
pling. We estimate a completeness limit of 10 Mq, cal- 
culated by determing the mass of a 3-(t detection at the 
upper distance limit of 5 kpc, assuming the dust temper- 
ature of the detection equals that of the sample average 
dust temperature, with T'd„st=44 K. The best fit power- 
law to the mass spectrum is fairly flat below ~ 80 Mq, 
with A''(> m) oc m~^'^'^. A break in the spectrum is seen 
around 100 Mq, above which the mass spectrum steepens 
to A^(> m) (X m~^^^'^. The mass distribution breakpoint 
is found comfortably above the completeness limit, sug- 
gesting this is not an observational artefact. 

We compared the index of these power-law fits to other 
observations, finding the distribution below 100 Mq to be 
significantly flatter than that seen in other studies. This 
points to either an absence or accelerated evolution of 
the lower-mass clumps, both of which would reduce the 
submm emission observed in the lower mass range. Most 
likely, this reflects an absence of low-mass clumps as we 
have only observed high-mass candidates, preferentially 
sampling only the high-mass tail of the initial mass dis- 
tribution, and additionally many close low-mass clumps 
undoubtedly lie unresolved. Correcting for these factors 
could raise the power-law index to a more typical IMF- 
like power-law slope, whereas at the moment the most 
similar power-law index for this region is that of clumps 
within molecular clouds, where the index averages around 
-0.60 (Kramer et al. 1998; Kramer et al. 1996). This is 
still a factor of four greater than our study, and without 
knowing how far we can minimize this difference we must 
question their true degree of similarity. 

Above the 100 Mq breakpoint, the slope of our sample 
becomes very similar to that of the fleld star IMF {N{> 
m) oc m-i-^ for tti/Mq > 0.5; Kroupa 2001). If the IMF- 
like distribution is valid, the apparently similar power- 
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Fig. 12. A histogram of the object mass, assuming the near kinematic distance (left-hand and centre plot) and far kinematic 

distance (right-hand plot) for distance- unresolved IRAS fields. The contribution of distance resolved detections to each histogram 
is plotted by an outlined histogram. The mass of each detection is calculated from 850 /xm continuum emission using a 100:1 gas- 
to-dust ratio. The centre panel display an expanded view of the near kinematic distance projection for clumps with M < SOOMq . 
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Mass Range 


Contribution (%) 


Type 


Mq 


Population 


Mass 


Brown dwarfs 


0.01-0.08 


37 


4.3 


M dwarfs 


0.08-0.5 


48 


28 


K dwarfs 


0.5-1.0 


8.9 


17 


Intermediate 


1.0-8.0 


5.7 


34 


stars 


> 8 


0.37 


17 



incomplete sampling 



Table 5. Stellar IMF, taken from Kroupa (2001). The third 
and fourth columns tabulate the contribution of each object 
type to the population and mass total. 
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Fig. 13. Cumulative mass distribution of all detections found 
at a kinematic distance of < 5 kpc, incorporating 65 clumps 
when assuming the near kinematic distance (upper curve) and 
32 clumps for the far kinematic distance (lower curve). The 
mass of each object is calculated from the 850 /jm emission 
using a 100:1 gas-to-dust ratio. The error bars correspond to 
Vn counting statistics. The thick dotted and thick dashed lines 
plot the lines of best fit for the first seven bins and the last four 
bins respectively. 



law index of main-sequence and massive pre-stellar clumps 
would suggest the star formation efSciency within these 
protoclusters is relatively mass invariant; it would appear 
that just the breakpoint shifts to lower mass as the core 
fragments and additional stars form. 

Assuming the multiple power-law IMF of Kroupa 2001, 
we may estimate the number of sources within an average 
350 Mq clump. Even if the star-formation efficiency e for a 
IOM0 star is as low as 5%, this still leaves sufficient mass 
to form an additional ^ 100 lower-mass objects created 
with efficiency e = 30%, distributed via the IMF number 
ratios given in Table 5. With a typical clump radius of 
~ 0.5', this would result in a typical stellar volume density 
of ^ 150 stars pc~^. 

With a median mass of ~143 Mq, the typical median 
clump could easily form at least one > 8 Mq star assum- 
ing an average 50% star-formation efficiency and an IMF- 



like mass distribution, assuming a high-mass star accounts 
for 17% of the total cluster mass (Table 5). However, 
the limited mass reservoir of the lowest-mass companion- 
less clumps means the star-formation efficiency e within 
these objects must be fairly high, as the total mass of 
gas and dust is close to that of a high-luminosity proto- 
star. The limited residual mass reservoir would also mean 
that not many lower-mass protostars can co-exist within 
the clumps. As a result, the stellar mass spectrum within 
these clumps must be skewed compared to the IMF of 
field stars or that within higher mass clumps. Although 
an isolated high-mass protostar has yet to be found, we 
suggest these clumps form ideal candidates of high-mass 
protostars with a minimal number of stellar companions. 

6. Discussion 

6.1. Comparison to 1.2mm observations 

The 850 /xm and 450 /xm emission we observe with the 
JCMT is found to be similar to the 1.2 mm IRAM contin- 
uum observations of Beuther et al. (2002a). The detections 
are comparable in size and morphology, although the ex- 
tended IRAM field of view reveals more sub-clumps per 
target field, strengthening the observation that high- mass 
star-forming clumps do not exist in isolation. 

The mass of coincident 850 /xm detections and 1.2 mm 
detections correlate well when the same gas and dust char- 
acteristics are assumed, as can be seen in the upper plot of 
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Fig. 14. A plot displaying the correlation between mass (upper 
plot) and peak flux (lower plot) of coincident IRAM 1.2mm 
detections and JCMT 850 jim detections. The dashed line in 
the lower plot traces the 1:1 mass ratio. 



Figure 14. Points lying away from the main trend gener- 
ally represent sources with companions lying in extended 
emission, suggesting the difference originates in the dif- 
ferent techniques used to measure the integrated emis- 
sion around extended and multiple component sources. 
Admittedly, the choice of where once source ends and an- 
other begins can be subjective, and we do not consider 
this a cause for concern. 

There is good agreement in the peak flux of 850 /um 
and 1.2 mm detections, which implies there is no large 
optical depth gradient between these wavelengths, and 
additionally that the studies are calibrated consistently 
with respect to one another, measuring approximately 7 
Jy/14.4" beam at 850 /im per 1 Jy/11" beam measured 
at 1.2 mm (lower plot of Figure 14). This is another in- 
dication that the 850 /xm and 1.2mm observations are de- 
tecting the same material and thus trace the same amount 
of mass. By confirming the mass of the clumps, we may 
place more confidence in the mass-luminosity relationship 
derived by Sridharan et al. (2002). 

6.2. Distance ambiguities 

We remain conscious that the distance ambiguities re- 
maining towards a number of sources could alter the con- 



clusions we are able to draw. Attempts to resolve the un- 
certainties using the C^^S(2-1) and CS(2-1) linewidths of 
Bcuther et al. (2002a), comparing the continuum-derived 
mass estimate with the virial mass at both kinematic dis- 
tances, proved fruitless. However, the simulations of Walsh 
et al. (2001) require the majority of our distance unre- 
solved sources to be placed at the near kinematic distance 
if we are to obtain M/L ratios similiar to those expected 
for the majority of 10^ — 10"^ M© clusters. 

The greybody analysis of SBSMW also gives an esti- 
mate of the bolometric luminosity of each IRAS source. 
Assuming oc L'^'^i^ (where Af* is the mass of the pro- 
tostar and L^oi is the SBSMW luminosity), similar to 
the mass-luminosity relationship for massive stars on the 
main sequence, it is simple to calculate the mass of the 
illuminating source. In comparing to the mass of the 
clump at both the near and far kinematic distance projec- 
tions, we found the luminosity of detections #64 (IRAS 
18553-F0414) and #74 (IRAS 19266-M745) to be incom- 
patible with their near-distance clump mass, suggesting 
that these sources are actually located at the far kine- 
matic distance. 

7. Conclusions 

We observed a sample of candidate high-mass protostars 
with the JCMT. Dust continuum emission was detected 
towards all sources, and from analysis of the 850 /um and 
450 /um maps we reach the following conclusions: 

— The average clump has a mass of ~ 330 M©, a total 
hydrogen column density of ~ 9 x 10^^ cm^^, and an 
angular diameter of 30", which projected to the mean 
projected distance of 4 kpc equals a linear diameter 
of 0.6 pc. Assuming spherical symmetry, these values 
translate to an average hydrogen number density of 
3 X 10^ cm^'^ in a typical clump. 

— The clumps we detect vary in mass from around 1 M© 
to over 1000 Mq. The high luminosity and low mass 
of the smallest clumps suggests they will only form a 
minimal number of stellar companions, and that we 
are close to probing the final mass of the most lumi- 
nous protostar. These small, isolated clumps may be 
analogues of the low- mass NII3 cores (Myers & Benson 
1983), from which we have learned a great deal about 
low mass star formation. As such, these small high lu- 
minosity cores may be the most promising sites for fur- 
ther high-resolution observations. The extremely large 
mass of the largest clumps suggest they may be form- 
ing proto-groups or proto-clusters. 

— Above 100 M0, the mass spectrum of the submm de- 
tections displays a distribution very similar to that 
of the field star IMF, suggesting the clump-mass to 
stellar-mass transfer efficiency is relatively mass invari- 
ant within these protoclusters. 

— A large fraction of the submm clumps we detect are 
roughly coincident with IRAS and MSX detections in 
the region. This suggests the short wavelength photons 
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and submm photons have the same origin, despite the 
typical envelopes being optically thick at wavelengths 
below 90 /im. The coincident detections could be con- 
ceivably be reconciled by non-uniformity of the enve- 
lope, eg. where outflow cavities provide a low-opacity 
escape route for short wavelength photons. 

— We also find some submm detections are significantly 
offset from their associated MSX and IRAS counter- 
parts, in agreement with Sridharan ct al. (2002). This 
raises the question of whether these objects are re- 
lated, or whether they are actually unique objects with 
different submm and IR properties. If the clumps we 
detect are potential proto-groups or proto-clusters, it 
is possible that young stars in the less dense, more 
transparent extremeties of the envelope may be more 
visible at shorter wavelengths. This possibility could 
explain both coincident and non- coincident long/short 
wavelength detections. 

The degree of coincidence may point to evidence of 
evolution, from the oldest, most evolved sources with 
MSX detections and no coincident submm flux, to 
MSX detections with some degree of submm flux, to 
the youngest, most embedded sources with potentially 
no MSX detection and high IR optical depth. Further 
investigation of these possibilities will require high res- 
olution observations at submm and far infra-red wave- 
lengths. 

— We find a companion clump fraction of ~ 0.6, empha- 
sising that clumps bearing high-mass stars do not form 
in isolation, and that they may lie in close proximity 
to other potentially star-forming clumps. 

— We measure the mean surface density of companions, 
and find the clump spatial density distribution peaks 
at a separation of r ^ 0.4 pc. The projected distribu- 
tion {Nproj{r) oc r^^-^) corresponds to a volume den- 
sity distribution such that N^oi{r) oc r"" "^^. 

— The mean spectral index of the dust emission, a, at 
the position of peak submm emission is 2.6 ± 0.4, and 
we observe both a morphologies both correlated and 
anti-corrclatcd with submm intensity. 

Peaked a morphologies can be reproduced with simple 
internal heating of a dusty envelope, with the observed 
a gradient resulting from the corresponding tempera- 
ture gradient through the envelope. 
An a-dip morphology could occur if the inner region of 
the clumps are cool with respect to their surroundings 
or if substantial grain growth occurs in these dense 
central regions. Although we suggest grain growth is 
the most likely factor, identifying the dominant mech- 
anism will require high-resolution observations of tem- 
perature tracers towards the regions. 

— The mean spectral index of the optical depth, (3, at the 
position of peak submm emission is 0.9 ± 0.4. This is 
lower than (3 towards in most dusty regions, but could 
be consistent with grain growth occuring in the dense 
clumps we observe. 
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